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ABSTRACT 

We extend a previous cross-calibration study by the International Astronomical Consortium for High Energy Calibration (IACHEC) 
on XMM-Newton/EPIC, Chandra/ ACIS and BeppoSAX/MECS X-ray instruments with clusters of galaxies (Nevalainen et al., 2010) 
to Suzaku/XIS instruments. Our aim is to study the accuracy of the energy-dependent effective area calibration of the XIS instru- 
ments using observations during 2005-2008. By comparison of the spectroscopic temperatures, fluxes and fit residuals obtained with 
Suzaku/XIS and XMM-Newton/EPIC-pn for the same cluster, we evaluated the systematic uncertainties of the energy dependence 
and the normalisation of the effective area between different detectors. The temperatures measured in the hard 2.0-7.0 keV energy 
band with all instruments are consistent within ~5 %. However, temperatures obtained with the XIS instruments in the soft 0.5-2.0 
keV band disagree by 9-29 %. We investigated residuals in the XIS soft band, which showed that if XIS0 effective area shape is 
accurately calibrated, the effective areas of XIS1 and XIS3 are overestimated below ~1.0 keV (or vice versa). Adjustments to the 
modelling of the column density of the XIS contaminant in the 3-6 arcmin extraction region while forcing consistent emission mod- 
els in each instrument for a given cluster significantly improved the fits. The oxygen column density in XIS 1 and XIS3 contaminant 
must be increased by ~l-2 x 10 17 cirT 2 in comparison to the values implemented in the current calibration, while the column density 
of the XIS0 contaminant given by the analysis is consistent with the public calibration. XIS soft band temperatures obtained with 
the modification to the column density of the contaminant agree better with temperatures obtained with the EPIC-pn instrument of 
XMM-Newton, than with those derived using the Chandra- ACIS instrument. However, comparison of hard band fluxes obtained using 
Suzaku-XIS to fluxes obtained using the Chandra-ACIS and EPIC-pn instruments proved inconclusive. 

Key words, instrumentation: miscellaneous - techniques: spectroscopic - galaxies: clusters: intracluster medium - X-rays: galaxies: 
clusters 

ON 

• , 1. Introduction high energy calibration and to supervise cross-calibration be- 

T ~* tween different observatories. This paper is based on the activi- 

O ■ Clusters of galaxies are bright, stable and extended objects and ties of ±e IACHEC clusters working group . The working group 

CO ; thus suitable for X-ray calibration. The X-ray spectra of clusters performed a study of clusters of galaxies in order to investi . 
ryH , consist typically of a combination of thermal bremsstrahlung and gate the effectiye ^ cross . calibration status bet ween XMM- 
> ■ colhsionally excited line emission. Temperature measurements NewtQn f Epic chandra f Aas and BeppoS AX / M ECS instru- 
of clusters above kT ~ 2 keV are dominated by the shape of the ments (Nevalainen et aL 20 10). They found that while the energy 
X-ray continuum. Consequently, by comparing temperatures of dependence of me effective area is accurate l y calibrated in the 
b ; the same cluster measured with different instruments, one can hard band (2 Q _ J Q tey) ±e eflfective ^ normalisa . 
get an indication on the cross-calibration status of the energy de- tions disagree by 5 _ 10 % . They also found significant cross- 
pendence of the effective area (defined here as the product of calibration differ ences in the soft energy band (0.5 - 2.0 keV). 
the mirror effective area, detector quantum efficiency and filter Th determined that tem peratures measured with Chandra / 
transmission) between the instruments. Similarly the normah- Aas m systematically „ 2 % higher than those measured with 
sation of the effective area affects the determination emission X MM-Newton / EPIC instruments. This amounts to a ~ 10 % 
measure and cluster flux, and comparison of cluster fluxes mea- uncertaint of the ACIS and EPIC soft band cross-calibration of 
sured with different instruments thus yields information about the effectiye ^ depe ndence. D ue to the high statistical 
the relative normalisation of the effective areas of the instru- weight of the soft band data> Nevalaine ]^^i] (gOlO) estimated 
ments. Analysis of the residuals of the spectral fits yields more that chandra ; ACIS and XMM-Newton / EPIC cluster temper- 
detailed information on the cross-calibration. atures measured in a typical 5 _ 7 keV energy band suffer 
The International Astronomical Consortium for High Energy from a systematic uncertainty of 10 - 15 % due to the cross- 
Calibration (IACHECU was formed to provide standards for calibration uncertainties 
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We extend the IA CHEC cluster work to include Suzaku / XIS 
dKovama et alJl2007h . We will evaluate the cross-calibration ac- 
curacy of the XIS instruments by comparison of measured tem- 
peratures and spectral fit resi duals. By com paring the measure- 
ments with XIS and EPIC-pn ( Striider et al. 2001) we aim to im- 
prove the knowledge about the soft band cross-calibration uncer- 
tainty between EPIC and ACIS.0 

This paper is organised as follows. In Section 2 we present 
the cluster sample used for our study, and in Section 3 we dis- 
cuss the selected spectral extraction regions. We present data re- 
duction and background modelling methods in Sections 4 and 5, 
respectively. We discuss the spectral analysis method in Section 
6 and present the results from the analysis in Sections 7 and 8. 
In Section 9, we investigate if uncertainties in the calibration 
of the optical blocking filter (OBF) molecular contaminant may 
affect the cross-calibration of the effective area energy depen- 
dence. Finally, we present our concluding remarks in Section 
10. Throughout this paper, we use 68 % (lcr) confidence level 
for errors, unless stated otherwise. 

2. Cluster sample 

We constructed a cluster sam ple consists of HIFLUGS clusters 
dReiprich & Bohringer 2002) with XMM-Newton and Suzaku 
observations publicly available in June 201 1, using our selection 
criteria detailed below. The clusters passing our selection criteria 
(see below) have been observed during the years 2005 - 2008 
(Suzaku) and 2000 - 2009 (XMM-Newton). The chosen clusters 
are all bright (flux ~ 10~ 12 - 10~ 10 erg / cm 2 / s in a 2 - 7 keV 
band) and nearby (z < 0.0753). General properties of the cluster 
sample are presented in Table Q] and [2] 

Our sample contains only such pointings where the off-axis 
angle between the cluster centre and the FOV centre is smaller 
than 1 arcmin. This ensures that we are studying the same de- 
tector regions and thus do not fold in any additional calibration 
uncertainties that depend on the position of the detector. 

We performed spectral analysis in two energy bands: 0.5 - 
2.0 keV (soft band) and 2.0 - 7.0 keV (hard band). We required 
that the observations have a minimum of 5 000 and a maximum 
of 500 000 data counts in the soft and hard energy bands. The 
lower limit ensures sufficient number of counts to obtain reliable 
statistics, while the upper limit eliminates such high signal-to- 
noise data where other calibration uncertainties, like gain and 
energy redistribution, are significant. We further excluded clus- 
ters where the total sky and non X-ray background exceeded 10 
% of the cluster flux (see Section for details on background 
modelling). 

In order to minimise additional uncertainties in the soft en- 
ergy band, we excluded clusters with a galactic absorption col- 
umn density Nh > 6 x 10 20 cirT 2 from the soft band sample (see 
Table [3}. Thus the soft band sample contains a smaller number 
of objects than the hard band sample. 

3. Extraction regions 

As we attempted to minimise PSF scatter and background con- 
tamination, we chose to use 3-6 arcmin annuli around the clus- 
ter centres as extraction regions (see Fig. [T}. By excluding the 
inner 3 arcmin we minimise the Suzaku PSF scatter from the 
cool core (see Section 15.5b . The width of the annulus is larger 
than the Suzaku PSF so that the PSF scatter to and from the 

2 Due to the limited FOV of ACIS-S we cannot perform a one-to-one 
comparison between the ACIS measurements and XIS. 



extraction region is minimised. By excluding the regions out- 
side 6 arcmin distance from the cluster centre, i.e. maintaining in 
the bright central regions, we minimize the effect of background 
contamination. 

The useful detector area (when excluding the CCD gaps and 
the regions around the calibration sources) varies somewhat be- 
tween different XIS instruments. The fraction of the full extrac- 
tion annuli covered by the useful detector area varies at most by 
12 % between the different instruments. We thus assume that the 
temperature measurements are not significantly affected due to 
the above difference in the covered regions. 

4. Data reduction 

The general properties of the observations are listed in Table|2] 

4.1. Suzaku 

We used HEASoffj release 6.11, containing Suzaku ftools ver- 
sion 18, for the processing. XIS calibration information was 
from XIS CALDB release 201 106080 and XRT calibration in- 
formation from XRT CALDB release 20080709.0 

We used the aepipeline tool to reprocess and screen the data 
and the xselect tool to extract spectra, images and light curves. 
We used standard filtering, with the exception that we eliminated 
times with geomagnetic cut-off rigidity (COR2) less than 6 and 
time since the South Atlantic Anomaly (SAA) passage less than 
436 sec. We excluded regions contaminated by the 55 Fe cali- 
bration sources and circular regions with a radius of 1 arcmin 
centred on bright point sources seen in XMM-Newton images 
(see Fig.Q}. 

We constructed energy redistribution and auxiliary response 
files with the xisrmfgen and xissimarfgen tools respectively. For 
each cluster, we used the XMM-Newton / EPIC-pn count rate 
map to weight the position-dependent effective area with the 
brightness distribution when producing the auxiliary responses. 
For the sky background modelling we produced the auxiliary re- 
sponses assuming a circular source of uniform brightness with a 
radius of 20 arcmin. We used makepi files version 20070124 and 
20100929, RMF parameter file version 20080901 and contami- 
nation file version 20091201 as included in XIS CALDB release 
20110608. 

We checked XIS 1 gain by searching for residuals in the spec- 
tral fits (Section [6) around the Fexxv lines. We found evidence 
for residuals in the Centaurus cluster and performed a gain fit in 
Xspec. However, the cluster temperature was unaffected. 

4.2. XMM-Newton EPIC-pn 

We used SAS version xmmsas_20110223_1801-l 1.0.0 and the 
latest calibration information in June 201 1 for EPIC-pn data re- 
duction. We used the epchain tool with the default parameters 
to produce event files. We filtered the event files, excluding bad 
pixels and CCD gaps (i.e. using expression flag=0), including 
only patterns 0-4. We generated simulated out-of-time event 
files, which we used to subtract events registered during CCD 
readout according to pn mode (Table|2]i. 



3 http://heasarc.nasa.gov/lheasoft/ 

4 http://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/data/suzaku/xis/ 
index/cif _suzaku_xis_20 1 1 0608 .html 

5 http://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/data/suzaku/xrt/ 
index/cif_suzaku_xrt_20080709.html 
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Table 1. General information on the cluster sample 



Name 


N H fl 


z" 


RA (J2000) c 


Dec (J2000) r 




[10 20 cm" 2 ] 




hh mm ss 




A 1060 


4.90 


0.0126 


10 36 42.6 


-27 31 44.2 


A 1795 


1.19 


0.0625 


13 48 53.0 


26 35 25.0 


A262 


5.67 


0.0163 


01 52 46.0 


36 09 09.1 


A3 112 


1.33 


0.0753 


03 17 57.7 


-44 14 18.3 


A496 


3.76 


0.0329 


04 33 38.1 


-13 15 40.0 


AWM7 


8.89 


0.0172 


02 54 27.4 


41 34 46.8 


Centaurus 


8.30 


0.0114 


12 48 48.9 


-41 18 46.9 


Coma 


0.87 


0.0231 


12 59 35.7 


27 57 34.0 


Ophiuchus 


19.3 


0.0280 


17 12 27.9 


-23 22 08.2 


Triangulum 


11.5 


0.0510 


16:38:18.6 


-64:20:37.8 



Notes. (o) LAB weighted average (Kalberla et al. 2005). {b) From NASA Extragalactic Database (http://ned.ipac.caltech.edu/! (c) Adopted coordi- 
nates of the cluster centre. 



Table 2. General information on the observations of the sample 







XMM-Newton 






Suzaku 








Name 


ID 


Obs.date 


Exposure time (ks) ° 


pn mode b 


ID 


Obs.date 


Exposure time (ks) " 






yyyy-mm-dd 


pn 






yyyy-mm-dd 


XIS0 


XIS1 


XIS 3 


A 1060 


0206230101 


2004-06-29 


14 


F 


800003010 


2005-11-22 


30 


30 


30 


A 1795 


0097820101 


2000-06-26 


34 


F 


800012010 


2005-12-10 


10 


10 


10 


A262 


0109980101 


2001-01-16 


16 


E 


802001010 


2007-08-17 


33 


33 


33 


A3 112 


0603050101 


2009-07-21 


63 


E 


803054010 


2008-05-23 


52 


52 


52 


A496 


0506260401 


2008-02-18 


33 


F 


803073010 


2008-08-02 


31 


31 


31 


AWM7 


0605540101 


2009-18-16 


76 


E 


801035010 


2006-08-07 


17 


17 


17 


Centaurus 


0406200101 


2006-07-24 


52 


E 


800014010 


2005-12-27 


30 


30 


30 


Coma 


0153750101 


2001-12-04 


17 


F 


801097010 


2006-05-31 


166 


166 


166 


Ophiuchus 


0505150101 


2007-09-02 


10 


E 


802046010 


2007-09-24 


75 


75 


75 


Triangulum 


0093620101 


2001-02-17 


3 


E 


803028010 


2008-10-11 


64 


64 


64 



Notes. {a) Clean exposure times after filtering, (A) F = Full Frame, E = Extended Full Frame. 



We used the evselect tool to extract spectra, images and light 
curves from the event file, excluding regions contaminated by 
bright point sources. We used the light curves to further filter 
the events by excluding the periods when the E > 10 keV band 
flux exceeded the quiescent rate by ~ 20% We used the rmfgen 
and arfgen tools to produce energy redistribution and auxiliary 
response files respectively. We used the arfgen tool in extended 
source configuration and weighted the response with an XMM- 
Newton image of the cluster in detector coordinates, binned in 
0.5 arcmin pixels. 

5. Background modelling 

As our sample consists of nearby clusters that fill the whole FOV, 
local estimates for the background emission are not available. 

5. 1 . Sky background 

We used the ROSAT All Sky Survey based HEASARC X-ray 
background tool0 to estimate the spectra of the sky background 
for both XMM-Newton / EPIC and Suzaku / XIS instruments. 
We extracted spectra of annular fields centred on the clusters, 
with an inner radius of 0.5 degrees and outer radius of 1 .0 de- 
grees. We fitted the spectra with a model consisting of an ab- 
sorbed and unabsorbed thermal MEKAL emission component 
with solar abundances (Mew e et al. I fl99l . describing emission 
from the Galactic halo and Local Hot Bubble respectively, and 

6 http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/xraybg/xraybg.pl 



an absorbed power-law with the index fixed to 1.45, describ- 
ing the extragalactic X-ray background^ We fixed the absorp- 
tion column densities to the value of the respective clusters (see 
Section [6). We added the local sky background emission model 
to all subsequent spectral fits. In order to account for the uncer- 
tainty in the sky background model, we allowed the temperatures 
and normalisations of the emission components to vary within 
l.Ocr of their best-fit values during subsequent spectral fits. 

5.2. XMM-Newton non X-ray background 

We extracted and investigated full FOV hard energy band (10 - 
14 keV) EPIC pn light curves to determine the quiescent non-X- 
ray background levels and performed GTI filtering by excluding 
flares. We modelled the quiescent par ticle background by sc al- 
ing a sample of closed filter spectra dNevalainen et al.ll2005l) to 
match hard band full FOV spectra. We removed the scaled closed 
filter spectra from the cluster spectra by using them as correction 
files is Xspec in subsequent pn spectral fits. 

We studied the uncertainty of the adopted particle back- 
ground model by varying the scaling of the closed filter spectra 
by + 10 % and interpreted the change in best fit temperature as 
a systematic effect due to background uncertainties. This had on 
average a ~ 3 % effect on the best fit temperatures in the hard 
band; in the soft band, the effect was less than 0.1% for all clus- 
ters in our sample. We propagated this uncertainty by adding it 
in quadrature to the statistical uncertainties of the temperatures. 

7 See e.g. http://heasarc.gsfc.nasa.gov/Tools/xraybg_help.html for 
references. 
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A1060 " "* A1795 




A A f — «. f~\ 50 100 150201 

49b 

Fig. 1. XIS1 raw count rate maps (not corrected for the background or vignetting) of the soft band cluster sample in a 0.5 - 7.0 keV 
energy band binned to 8 pixels (~ 0. 14 arcmin) per bin. The inner and outer radii of the 3-6 arcmin extraction region are denoted 
by white circles, the green circles show the excluded regions centred on point sources seen in XMM-Newton EPIC / pn images. 
Scaling of the maps are adjusted for plot clarity and do not reflect the relative brightness of the clusters. The colour bar shows the 
the total number of counts per bin. 



5.3. Suzaku non X-ray background 

We used the xisnxbgen tool with default settings to generate 
spectra of the Suzaku XIS non X-ray background (NXB) based 
on night Earth observations. We used COR2 weighting for NXB 
construction and to verify the normalisation of the NXB gener- 
ated using the xisnxbgen tool, we compared the > 10 keV NXB 
flux with that in the cluster data and found them consistent. We 



subsequently subtracted the NXB spectra from the cluster spec- 
tra in all Suzaku XIS spectral fits by using the NXB spectra as 
background files in Xspec and thus considered the uncertainty of 
the NXB in the spectral fits. 
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XIS1 A496 3 - 6 arcmin 



10.0000 



1.0000 



0.1000 



c 0.0100 



0.0010 



0.0001 




Energy [keV] 



Fig. 2. The total emission (blue crosses), sky background (black 
curve) and non X-ray background (red crosses) emission of 
A496 in the 3-6 arcmin region for XIS1. The sky back- 
ground is modelled from ROS AT All Sky Survey data and folded 
through the XIS1 response. The XIS1 observation of A496 has 
the strongest relative background emission of the cluster sample. 



5.4. Charge exchange contamination 

Geocoronal and heliosphe ric Solar wind charge e xc hange 
(SWCXt (IWargelin et all l2004t ISnowden et aD l2004t 
iFuiimoto et al.l I2007I) may produce contamination to the 
cluster emission in terms of emission lines in the 0.5 - 0.9 
keV band. The total flux of the SWCX varies by an order 
of magnitude at a timescale of hours or minutes, reaching 
a maximum level comparable to that of the sky background 
(Galactic emission + CXB). To obtain an conservative estimate 
of the possible SWCX effect on our results, we examined the 
variation of the best-fit 0.5 - 2.0 keV energy band temperatures 
of XIS 1 when keeping the sky background model fixed to RASS 
values, and by increasing the normalisation by 100 %. Due to 
the low level of the sky background compared to the cluster flux 
in our sample the above effect was insignificant. 

5.5. XIS point spread function 

Cluster flux originating from a given extraction region may be 
conta minated by PSF scattered flux from neighbouring regions 
(e.g. ISato et ai1l2007l) . Since Suzaku XIS has a relatively large 
PSF (~ 2 arcmin half-power diameter), scatter from cool cores 
may affect observed XIS temperatures. Thus we examined how 
well PSF scatter to and from our extraction region (3-6 arcmin 
annulus, see Section [3} is minimised by the usage of extraction 
regions wider and further away from the centre than the XIS PSF. 

The exact PSF scattering effects depend on the brightness 
of the cool centre, the shape of the surface density profile, the 
temperature of the cluster and the cluster redshift, and therefore 
they can vary from cluster to cluster. We choose to use A 1795 
and A3 112 as conservative examples for the PSF-induced sys- 
tematic uncertainty, as they possess cool cores which are an or- 
der of magnitude b righter than the cores of the rest of the sample 
(Pere s et al 1fl998h . 

We performed ray-tracing simulations of A 1795 with xis- 
sim to estimate the fraction of phot ons that originate from the 
central cool region (r ~ 1.5 arcmin; Vikhl inin et al.l (2005)) but 
are scattered into the 3-6 arcmin annulus. We used an expo- 



sure corrected 0.5 - 7.5 keV energy band XMM-Newton EPIC 
/ pn rate map of A 1795 as a source brightness model and, since 
the energy depen dence of the PSF is known to be small (e.g. 
Sato et al] (120071) Appendix 1), performed the simulation at 1 
keV. We accounted for typical pointing errors by using the at- 
titude and exposure information from the Suzaku observations 
of A 1795 used in this work. The simulations show that about 
16 % of the flux detected in the 3-6 arcmin extraction region 
originates from the inner 1.5 arcmin cool core. There is a small 
variation of scattered light between detectors, with values of 13 
% (XIS0), 18 % (XIS1), and 15 % (XIS0). 

We then simulated a XIS 1 spectrum in the 3-6 arcmin an- 
nulus containing one MEKAL component for the emission orig- 
inating from that region and another for the PSF scattered emis- 
sion fro m the cool core. We u sed the Chandra measurements of 
A1795 (I Vikhlinin et al.l20 05) for the temperature and abundance 
of the two components: kT = 6.5 keV, 0.25 solar and kT = 4.75 
keV, 0.55 solar, respectively. We used relative model normali- 
sations of 0.84 and 0.16 respectively, given by the above PSF 
simulations. We fitted the simulated spectrum with a single tem- 
perature MEKAL model in the soft 0.5 - 2.0 keV and hard 2.0 - 
7.0 keV energy bands. The resulting best-fit temperatures were 
lower by ~ 7 % and ~ 6 % in the soft and hard energy bands 
respectively, compared to those obtained from the fits to the ob- 
ser vational data . 

iLehto et all (1201 Ol) performed an analogous XIS PSF scatter 
ray-tracing simulation on A3 112 data using a similar 3 - 6 ar- 
cmin extraction region as used for our study. They found that 
~ 10 % of the flux in the extraction region originates from the 
core, but best-fit temperatures were consistent within 2 % in all 
energy bands studied. The other clusters are weaker (or non-) 
cool cores, so we therefore expect lower biases for all of them. 
Thus we conclude that cool core contributions may bias Suzaku 
XIS temperatures low, but only up to 6 -7 % in the worst case, 
and much lower in most cases. However, as the three XIS in- 
struments have similar scattered light properties, they would be 
biased by similar amounts and the relative XIS temperature dif- 
ferences would be unaffected. Consequently, the effect is only 
limited to XIS to pn temperature comparisons. 

6. Methods 



iNevalainen et a l. (2010) found that several independent mis- 
sions / instruments (XMM-Newton / EPIC, Chandra / ACIS, 
BeppoSAX / MECS) and two independent physical processes 
(bremsstrahlung and collisionally excited line emission) yielded 
consistent temperature s in the 2-7 keV band. However, 
Nevalainen et al] (1201 Oh also found that the 2 - 7 keV band EPIC 
/ ACIS fluxes disagree by 5 - 10 %. Here we test the above find- 
ings with Suzaku / XIS by studying the cross-calibration of the 
the effective area by spectroscopic analysis of cluster data ob- 
tained with Suzaku / XIS instruments and the XMM-Newton / 
EPIC-pn instrument. 

6. 1 . Spectral analysis 

We extracted spectra of the clusters in our sample from annu- 
lar 3 - 6 arcmin regions centred on the clusters (see Section 
0. We binned the spectra to a minimum of 100 counts per 
bin. We fitted the NXB-subtracted spectra (see Sections 15.31 
and 15.2b with a model c onsisting of an ab sorbed single tem- 
perature MEKAL model (iMewe et alJll995|) and the sky back- 
ground c omponent (see Section [ 5.11 . We used metal abun- 
dances o f G revesse &Sauvall (ll998). We modelled the Galactic 
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absorption with the PHABS model, fixing the column den- 
sities to the values obtained with 21 cm radio observations 
(Kalberla et al. 200 5b and using the absorption cross sections of 
Balu cinska-Church & McCammonl(ll992l) . We allowed the tem- 
perature, metal abundance and normalisation of the cluster emis- 
sion model to vary in the fits. We performed the fits in two energy 
bands: 0.5 - 2.0 keV (soft band) and 2.0 - 7.0 keV (hard band). 

In addition to excluding clusters with Nh > 6 x 10 20 cm 2 
from the soft band sample (Section we excluded Coma from 
the soft band sample as we were not able to produce single tem- 
perature spectral fits of satisfactory statistical quality. Thus the 
soft band sample contains a smaller number of objects than the 
hard band sample. The properties of the observations are pre- 
sented in Table |3] 

We experimented with using C statistics (ICashll 19791) in ad- 
dition to standard x 2 statistics to estimate the best-fit parame- 
ters. We compared the temperatures obtained by fitting unbinned 
XIS1 spectra using C statistics to XIS1 temperatures obtained 
using x 1 statistics and binning the spectra to a minimum of 100 
counts per bin. Using the statistical analysis described below, we 
found that the C statistics yields systematically higher tempera- 
tures, with an average temperature difference of 0.5 % in the soft 
0.5 - 2.0 keV energy band (statistical significance 0.2cr) and 0.8 
% in the hard 2.0 - 7.0 keV energy band (0.8<x). Thus the choice 
of statistics in spectral fits can have an impact on cluster tem- 
peratures, but it has a low statistical significance and the effect 
is consequently minimal due to the high number of counts in 
our spectra. We will use binned spectra and x 1 statistics for the 
remainder of this paper. 

6.2. Comparison of temperatures and fluxes 

We computed the difference in temperature or flux measure- 
ments and their statistical significance for each cluster and in- 
strument pairing in terms of a fraction of the average value, 
f T = AT/ < T >= (r 1 -r 2 )/{(r I +r 2 )/2}and ( r r = cr AT / < T >; 
Sf = AF/ < F >= (Fi - F 2 )l{{Fi + F 2 )/2] and o> = cr AF / < 
F > respectively. As a systematic uncertainty in calibration 
would tend to drive the mean of the temperature or flux differ- 
ences away from zero, we attributed fj or Sf as a measure of the 
relative systematic uncertainty in the cross-calibration of a given 
instrument pair. 

6.3. Stacked residuals 

To examine the cross-calibration accuracy of the effective area of 
the studied instruments , we applied the stacked residuals method 
dLonginotti et al.ll2008l) . This method involves a choice of a ref- 
erence instrument, against which the other instruments are com- 
pared. Since we do not know which, if any, instrument has its 
effective area perfectly calibrated, the choice of a reference in- 
strument is arbitrary and the comparison is only relative. 

The stacked residuals method for a pair of instruments, II 
and 12, involves the following steps: 

- We start with the data (data nja) an d responses (resp 11,12), 
defined here as the energy response matrix (RMF) multiplied 
with the effective area (ARF), of the instruments studied. For 
this example we use II as the reference instrument and thus 
perform spectral analysis as explained above for each cluster 
to fit a reference model (model n). 

- In order to be able to compare the data of different instru- 
ments we regroup data 0,12 to equal binning 



- For each source, we convolve the reference model with the 
response of the instruments in order to produce an II -based 
model prediction (pred 0,12 = model n ® resp 11,12) 

- We divide the data of each cluster by the above model 
prediction to obtain residuals over the Il-based prediction 
(resid U j2 = data n ,i2 / pred i U2 ) 

- We divide the above residuals of 12 with II residuals, in order 
to remove effects of possible calibration uncertainties of II, 
obtaining ratio R21 for 12: 

R21 = resid xi I resid n 

data 12 model n <g> resp ,, 

— x (y\ 

model 11 ® resp I2 data n 

R21 is a useful measure of the effective area cross-calibration, 
since the value for each energy bin would be unity if the 
cross-calibration of II and 12 effective areas is consistent. 

- We calculate median and median absolute deviation of R21 
values for each energy channel for the cluster sample. 

7. Results 

The best-fit single temperature models, data and spectral param- 
eters are shown in Appendix lAl 

7.1. Hard band 

We found that the hard band (2.0 - 7.0 keV) temperatures mea- 
sured with XIS0 and XIS3 agree, with an average difference of ~ 
1 % (statistical significance 0.7cr). Temperatures measured with 
XIS1 are on average ~ 5 % (5.6<x) and ~ 5 % (4.9<x) higher than 
those measured with XIS0 and XIS3 respectively (see Fig. [3). 

The simultaneous Suzaku / XMM-Newton o bservations of 
BL Lac object PKS2 155-304 llshida et alJl20Tlh are useful for 
comparison with our cluster results since the observations cover 
the sam e time period (2005 - 2008). The XIS CALDB version 
used bv llshida et aH (l201ll) for PKS2155-304 (dated July 30, 
2010) differs from the one used for our cluster study (dated June 

8, 20 1 1 ) in the model used for the molecular OBF contamination 
(see Section 8). Thus there should be little difference in the cali- 
bration in the hard band. The best-fit flux ratios of PKS2 155-304 
indicate that XIS 1 yields somewhat harder spectra in the 2-7 
keV band than XIS0, consistent with our results. 

The hard band best-fit residuals of the clusters using the XIS 
instruments are rather consistent with unity (see Fig. |4} except 
that there are some systematic residuals around 3.0 - 3.5 keV 
energies at a few % level, especially in the XIS 1 data. These fea- 
tures are probably causing the small differences in XIS hard band 
temperatures (see above). This behaviour is possibly connected 
to the PSF effect, see below. 

In order to derive the effects of possible cross-calibration un- 
certainties in the shape of the effective area in the hard band be- 
tween Suzaku-XIS and EPIC-pn instruments, we compared the 
XIS hard band temperatures of our sample to those obtained with 
EPIC-pn. We found that the XIS1 and EPIC-pn temperatures 
are consistent within the statistical uncertainties: on average the 
XIS 1 temperatures are ~ 2 % ( 1 .7<x) lower than the EPIC-pn val- 
ues. XIS0 and XIS3 temperatures are respectively on average ~ 
6 % (4.9<x) and ~ 5 % (4.2cr) lower than pn temperatures (see 
Fig. [5]), indicating some cross-calibration problems. However, 
the lower XIS temperatures may be due to the uncorrected PSF 
scattering from the cool core which may bias XIS temperatures 
low by a maximum of ~ 6 % in the hard band (see Section l53T l. 
Indeed, the stacked XIS / pn residuals indicate some excess data 
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Table 3. Properties of the observations 



Name 


Background" / total flux 




Total data counts 






pn 


XIS0 


XIS1 


XIS3 


pn 


XIS0 


XIS1 


XIS3 








Soft energy band (0.5 - 


2.0 keV) 






A 1060 


0.04 


0.03 


0.03 


0.03 


73473 


24723 


38605 


22610 


A1795 


0.04 


0.05 


0.06 


0.06 


149804 


8296 


13926 


7537 


A262 


0.05 


0.04 


0.05 


0.04 


40870 


14479 


21265 


13668 


A3112 


0.08 


0.07 


0.07 


0.08 


95084 


10865 


17928 


11140 


A496 


0.02 


0.08 


0.10 


0.09 


134661 


20235 


29516 


18686 








Hard energy band (2.0 - 


- 7.0 keV) 






A 1060 


0.04 


0.03 


0.04 


0.03 


18789 


12033 


13982 


11898 


A1795 


0.02 


0.04 


0.04 


0.04 


48455 


5034 


6125 


5148 


A262 


0.07 


0.06 


0.06 


0.07 


8151 


5958 


6829 


6829 


A3112 


0.08 


0.06 


0.08 


0.06 


28344 


7014 


8628 


7298 


A496 


0.02 


0.06 


0.10 


0.06 


42088 


13524 


14770 


12786 


AWM7 


0.03 


0.02 


0.02 


0.02 


176775 


12786 


12235 


11008 


Centaurus 


0.05 


0.04 


0.06 


0.05 


115211 


18242 


20907 


17953 


Coma 


0.02 


0.02 


0.02 


0.02 


74184 


220176 


242263 


229227 


Ophiuchus 


0.01 


0.01 


0.01 


0.01 


114581 


265097 


287684 


271563 


Triangulum 


0.02 


0.01 


0.01 


0.02 


11799 


70564 


74603 


68502 



Notes. {a> total sky and non X-ray background 




XIS1 hord bond T [keV] XIS3 hord bond T [keV] XIS3 hord bond T [keV] 

Fig. 3. Upper panel: The best-fit hard band temperatures (asterisks) and 1 <x uncertainties (boxes) for Suzaku XIS instruments 
using the public calibration. The solid black line is an identity line, drawn as a reference. Lower panel: The relative temperature 
differences fj (diamonds) between different Suzaku XIS instruments, and their lcr uncertainties. The dotted and dashed lines show 
the weighted means of the relative temperature difference and corresponding uncertainties. 



at the lowest energies, reaching a maximum of ~10% at ~ 3 keV 
energies, which is possibly due to the PSF scatter from the cool 
core (see Fig. Since the effective area energy dependence of 
XMM-Ne wton' s EPIC-pn is most likely well calibrated in the 
hard band (Nevalai nen et al.ll2010h . the shape of the hard band 
effective areas of the XIS instruments are correctly calibrated 
within the statistical and PSF-related uncertainties, amounting 
to ~ 5% uncertainties on the measurement of the hard band tem- 
perature. 



7.2. Soft band 



The temperature measurement of hot clusters is not optimally 
achieved in the 0.5 - 2.0 keV band as the temperature depen- 
dence of the continuum shape is relatively weak since the cut-off 
is at higher energies. Due to weak line emission in the soft en- 
ergy band at typical cluster temperatures, possible uncertainties 
in metal abundance measurements mostly produce variation in 
continuum emission, which in turn may affect the temperature 
measurements. However, our aim here is to yield a phenomeno- 
logical, rather than physical model of the combination of cluster 
continuum emission and effective area calibration uncertainties. 
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Fig. 5. Upper panel: The best-fit hard band temperatures (asterisks) and 1 cr uncertainties (boxes) for XMM-Newton EPIC-pn and 
Suzaku XIS instruments using the public calibration. The solid black line is an identity line, drawn as a reference. Lower panel: The 
relative temperature differences f? (diamonds) between different XMM-Newton EPIC-pn and Suzaku XIS instruments, and their 
lcr uncertainties. The dotted and dashed lines show the weighted means of the relative temperature difference and corresponding 
uncertainties. 
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Fig. 4. The median and the median absolute deviation of the hard 
band residuals (data / best-fit model) for XISO (top), XIS 1 (mid- 
dle) and XIS3 (bottom). The model is an absorbed single tem- 
perature MEKAL whose temperatures and metal abundances are 
independent for each instrument, fitted to 2.0 - 7.0 keV band us- 
ing the public calibration. The data are regrouped to uniform 
binning of 0.5 keV. 
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Fig. 6. The median + median absolute deviation of hard energy 
band stacked residuals of XISO / pn Ro pn (top panel), XIS 1 / pn 
Ri pn (middle panel) and XIS3 / pn R.3 pn (bottom panel). The data 
are regrouped to uniform binning of 0.5 keV. 



Thus the soft band temperatures are still useful for the purposes 
of cross-calibration of the effective area. 

7.2.1. Temperatures and residuals 

We fitted the soft band spectra of the clusters in the soft band 
sample with an absorbed single-temperature MEKAL model 
(see Table lA.fl i. The soft band temperatures obtained using XIS 1 
and XIS3 differ by ~ 9 % , which is not very significant {2.9(f). 



Unfortunately XISO temperatures disagree significantly from the 
other two instruments. They are on average ~ 29 % (1 1.2cr) and 
~ 23 % (J Act) lower than XIS1 and XIS3 temperatures respec- 
tively (see Fig. I2]L 

The llshida etal] d201 lb work on PKS2 1 55-304 also indicates 
that XISO yields softer spectra than XIS1 in the 0.5 - 2.0 keV 
band during years 2005 - 2008, qualitatively consiste nt with our 
cluster results. However, as noted in Section 7.1, llshida etakl 
(1201 lh used an older XIS calibration with a different model for 
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Fig. 8. The median and the median absolute deviation of the 
residuals (data / best-fit model) of the soft band cluster sample 
for XISO (black), XIS1 (red) and XIS3 (blue). The model is an 
absorbed single temperature MEKAL whose temperatures and 
metal abundances are independent for each instrument, fitted to 
0.5 - 2.0 keV band using the public calibration. The data are 
regrouped to uniform binning of 127.75 keV. 
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Fig. 11. The median + median absolute deviation of the soft en- 
ergy band stacked residuals of XISO / pn R Upn (top panel), XIS 1 
/ pn Ri pn (middle panel) and XIS3 / pn R3 pn (bottom panel). The 
data are regrouped to uniform binning of 127.75 keV. 



the OBF contaminant, making a more rigorous comparison to 
that work difficult. 

XIS soft band residuals exhibit a systematic behaviour: 
they reach a maximum at 1.3 - 1.4 keV while they decrease 
with lower and higher energy (see Fig. |8). We attempted two- 
temperature modelling, but this did not significantly improve the 
spectral fits. Importantly, the residuals of the XISO unit, which 
yielded significantly lower soft band temperatures (see above) 
are rather consistent with unity below 1 keV energies, while 
those of XIS1 and XIS 3 decrease significantly below unity. 
Assuming for a moment that the shape of the XISO effective area 
below 1 keV is accurately calibrated, this indicates that the ef- 
fective areas of XIS1 and XIS3 are overestimated increasingly 
with lower energy. 

The comparison between XIS and XMM-Newton / pn soft 
band temperatures complicates the situation further: XISO values 
are on average ~ 14 % (6.2cr) lower than those of pn, whereas 
XIS 1 and XIS3 values are respectively on average ~ 18 % (7.4<x) 
and ~9 % (3. Oct) higher than those ofpn (see Fig . |9j . Since PSF 
scatter from the cool core could bias XIS temperatures low by 7 
% in only the worst case, this effect is unable to fully account for 
the XISO / pn temperature difference, and in fact any correction 
for it would exacerbate XIS 1 / pn and X IS3 / pn difference s. 

The observations of PKS2 155-304 dlshida et alj|201 lh also 
indicate that XIS1 yields harder spectra than pn in the 0.5 - 2.0 
keV energy band. As Chandra / ACIS soft band cluster temper- 
atures jtfej}n_average_~ 20% higher than pn soft band tempera- 
tures dNevalainen et al.l l201Q). XIS1 / ACIS are the only instru- 
ment pair that yields approximatively consistent sample average 
soft band temperatures. Due to these disagreements we do not 
have clear indications of which instruments are most accurately 
calibrated in the soft energy band. 

7.2.2. Stacked residuals 

However, this unfortunate situation does not prohibit us from 
deriving information about XIS cross-calibration by analysing 
the cluster measurements. The above XIS temperature and resid- 



ual discrepancies are indicative of uncertainties of the cross- 
calibration of the energy dependence of the effective area of 
different XIS instruments. As XIS1 as a back-illuminated (BI) 
device has the largest effective area of the three XIS CCDs, 
we adopted it as the reference instrument for a relative cross- 
calibration comparison between the XIS instruments using the 
stacked residuals method (see Section 16.3b . The analysis indi- 
cated a systematic difference between the cross-calibration of 
XISO / XIS 1 effective areas, which increases towards lower ener- 
gies (see Fig. [Toll. The median of XISO / XIS1 residuals reaches 
a value of ~ 1.2 at 0.5 keV, indicating ~ 20 % cross-calibration 
uncertainties. The case of XIS3 / XIS1 is less problematic, since 
the median of XIS3 / XIS1 differs significantly from unity only 
in the energy bin at ~ 0.5 - 0.7 keV. 

In order to evaluate the cross-calibration uncertainties of the 
soft band effective area shape between Suzaku / XIS and XMM- 
Newton / pn instruments, indicated by the differences of the soft 
band temperatures (see Section l7.2.H . we calculated the stacked 
residuals using EPIC-pn as the reference instrument (Fig. fTTb . 
While XISO / pn residuals are rather consistent with unity at the 
lowest energies, the data indicates that at energies above 1.3 keV 
XISO flux is lower than EPIC-pn flux by ~ 10 %. This effect 
would account for the systematically lower soft band XISO tem- 
peratures, when compared to EPIC-pn. On the other hand, XIS 1 
and XIS3 instruments yield ~ 10 - 20 % lower flux than EPIC-pn 
at energies below 1 keV. This results in higher soft band temper- 
atures for XIS 1 and XIS3, when compared to EPIC-pn. 

However, we do not have indications of which instrument is 
most accurately calibrated in the soft band. We thus conclude 
that if fitting only soft band data (0.5-2.0 keV), as is often done 
for high redshift clusters whose emission in the hard band is 
small compared to the background, using Suzaku-XIS or EPIC- 
pn instruments, the cross-calibration uncertainties between these 
instruments yield a maximum systematic uncertainty of ~ 30% 
on the derived temperatures. 
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Fig. 7. Upper panel: The best-fit soft band temperatures (asterisks) and 1 cr uncertainties (boxes) for Suzaku XIS instruments using 
the public calibration. The solid black line is an identity line, drawn as a reference. Lower panel: The relative temperature differences 
fj (diamonds) between different Suzaku XIS instruments, and their 1 cr uncertainties. The dotted and dashed lines show the weighted 
means of the relative temperature difference and corresponding uncertainties. 
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Fig. 9. Upper panel: The best-fit soft band temperatures (asterisks) and 1 cr uncertainties (boxes) for XMM-Newton EPIC-pn and 
Suzaku XIS instruments using the public calibration. The solid black line is an identity line, drawn as a reference. Lower panel: 
The relative temperature differences fj (diamonds) between XMM-Newton EPIC-pn and Suzaku XIS instruments, and their lcr 
uncertainties. The dotted and dashed lines show the weighted means of the relative temperature difference and corresponding 
uncertainties. 



7.3. Full band 

We then proceeded to evaluate the effects of the above cross- 
calibration uncertainties on standard astrophysical spectral anal- 
ysis of the cluster spectra, i.e. modelling and analysis of the full 
useful energy band of 0.5 - 7.0 keV, using the soft band sample. 
These results should be applicable for any type of source which 



produces such continuum-dominated X-ray emission in this en- 
ergy band which has significant statistical weight at energies be- 
low 2 keV. However, the usage of a singe-temperature emission 
model in the full band may not be physically justified due to 
possible multi-temperature structure of the intracluster medium. 
Indeed, while the XIS residuals show systematic features due to 
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Fig. 10. The median + median absolute deviation of soft energy band stacked residuals of XIS0 / XIS1 Roi (left panel) and XIS3 / 
XIS1 R31 (right panel) using the public calibration (dotted line) and the modified responses described in this work (solid line). The 
data are regrouped to uniform binning of 127.75 keV. 
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Fig. 12. The median and the median absolute deviation of the 
residuals (data / best-fit model) for the soft band cluster sample 
when fitting the full 0.5-7.0 keV band using the public calibra- 
tion for XIS0 (black), XIS1 (red) and XIS3 (blue). The model 
is an absorbed single temperature MEKAL whose temperatures 
and metal abundances are independent for each instrument. The 
data are regrouped to uniform binning of 0.5 keV. 



the calibration uncertainties discussed above, a part of the resid- 
uals may originate from multi-temperature structures (see Fig. 
[T2b . Yet the usage of the single-temperature emission model has 
the advantage of yielding a simple and useful measure for cross- 
calibration uncertainties. 



7.3.1. Temperatures 

While the cross-calibration uncertainties between XIS instru- 
ments and between Suzaku-XIS and EPIC-pn instruments were 
quite large below 2 keV (see Section 17721 . the largest discrepan- 
cies were concentrated in a rather narrow energy band of 0.5 - 
0.8 keV (see Figs. fTOlfTTI ). The emission of this band is relatively 
small when compared to the total emission of the full band and 



thus the effect of the soft band discrepancies on full band temper- 
atures is small. The full band analysis resulted in rather similar 
temperatures as the well calibrated hard bandfl 

Comparison of best-fit temperatures using the full band 
yielded that the effect of the lowest energy XIS cross-calibration 
uncertainties to the full band temperatures are ~ 3 %. The effect 
of Suzaku-XIS / EPIC-pn soft energy calibration uncertainties 
on full band temperatures is ~ 1 - 5 %, depending on the XIS 
instrument (see Fig.fTJb. 

7.3.2. Stacked residuals 

When using the stacked residuals method for evaluation of the 
cross-calibration uncertainties of the effective area, a very de- 
tailed physically sound model is not necessary: the usage of the 
stacked residuals formula (Eq.[T) will take into account and cor- 
rect for the deviations of the best-fit model and the data of the 
reference instrument. Assuming that the calibration of the ref- 
erence instrument is correct, the phenomenological model (i.e. 
the best-fit model corrected for the deviations with the data) de- 
rived for the reference instrument still yields a prediction which 
matches the data of the studied instrument, if the calibration of 
the studied instrument is accurate. 

Yet the stacked residuals for XIS instruments (using XIS 1 
as a reference as before) in the full band yield systematic dif- 
ferences at energies above 2 keV (Fig. fl4ll . inconsistent with 
the above argumentation and the fact that in the hard band we 
found a very good agreement between the different XIS instru- 
ments (Figs. [3] and |4j . We think that the reason for this is that 
the sample used for the hard band studies is larger than the one 
used for the soft and full band study due to our selection criteria 
(Section [2). In the hard band the dispersion of the residuals of 
individual clusters is bigger for the smaller sample where one or 
two deviant objects may cause the difference. Thus, the full band 
stacked residuals analysis above 2 keV energies should be taken 
with caution. 



8 The differences in the hard and wide band samples also play a role 
here, see Section l7.3.2l 
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Fig. 13. Upper panel: The best-fit full band temperatures (asterisks) and lcr uncertainties (boxes) for XMM-Newton EPIC-pn and 
Suzaku XIS instruments using the public calibration. The solid black line is an identity line, drawn as a reference. Lower panel: 
The relative temperature differences fj (diamonds) between XMM-Newton EPIC-pn and Suzaku XIS instruments, and their lcr 
uncertainties. The dotted and dashed lines show the weighted means of the relative temperature difference and corresponding 
uncertainties. 
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Fig. 14. The median + median absolute deviation of the full energy band stacked residuals of XISO / XIS 1 Roi (left panel) and XIS3 
/ XIS1 R31 (right panel) using the public calibration. Note that the cluster sample used here is different from the one used for the 
hard band. The data are regrouped to uniform binning of 0.5 ke V. 



A similar problem is present when analysing the stacked 
residuals of the XIS instruments using EPIC-pn as a reference in- 
strument: while the XIS / pn stacked hard band residuals showed 
a very good consistence (Figs. and |6), the stacked XIS / pn 
residuals of the smaller sample used for the full band (Fig. \T5[ 
exhibit systematic differences at energies above 2 keV. As above, 
we think that the differences of the samples are causing this dis- 
crepancy, and that the XIS / pn cross-calibration uncertainties 
above 2 keV are better evaluated using the larger hard band sam- 
ple. Thus, with the current data we cannot properly estimate the 



effect of the full band cross-calibration uncertainties between the 
Suzaku-XIS instruments and between the Suzaku-XIS and the 
EPIC-pn instruments. 

7.4. Suzaku old calibration 

In order to study the effects of the different calibration ver- 
sions on the discrepancies detected in the XIS data in both 
hard and soft energy bands, we reprocessed the data ad- 
ditionally with the previous Suzaku ftools version 17. The 
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Fig. 15. The median + median absolute deviation full band 
stacked residuals of XlO/pn (Ro pn , top panel), XlSl/pn (Ri pn , 
middle panel) and XIS3/pn (R3 pn , bottom panel). The data are 
regrouped to uniform binning of 0.5 keV. 



older version (XIS CALDB release 20110210) differs mainly 
by containing the previous model for the XIS contamina- 
tion (aejd*_contami^0081023.fits for XIS1 and XIS3, and 
ae_xi*_contami_20090128.fits for XISO) where the carbon-to- 
oxygen ratio was fixed to 6.0. 

The different calibration versions changed the best-fit hard 
energy band temperatures by less than 0.1% and the soft en- 
ergy band temperatures by less than 3.5 %. Thus, temperature 
discrepancies and consistencies reported above also apply to the 
old calibration. We will not study the effects of different XIS 
calibration versions in any more detail. 

8. Fluxes 

In addition to the cross-calibration of the shape of the effec- 
tive area, the cluster data can also be used to assess the cross- 
calibration of the normalisation of the effective area. To the first 
order the flux of the cluster emission depends linearly on the 
accuracy of the normalisation of the effective area. Due to the 
differences of the composition of the samples used for the anal- 
ysis of the hard and soft band (see Section|2]i, we will only study 
fluxes in the hard band here. 

We examined the hard band fluxes given by the best-fit un- 
absorbed MEKAL models presented in (Section lTTTI i. We calcu- 
lated the fluxes and their statistical uncertainties using the cflux 
model in Xspec. Due to masked point sources, bad columns and 
pn CCD gaps, the covered fraction of the 3-6 arcmin annulus 
used for spectral extraction varies between the instruments. To 
enable meaningful comparison of fluxes we calculated the flux 
in terms of mean surface brightness by dividing the total flux 
from the extraction region with the useful detector area for each 
cluster. Since the emission is not constant with radius, the dif- 
ference in covered areas might introduce some differences to the 
flux of a given cluster measured with different instruments. 

We found that the hard band fluxes measured with XIS in- 
struments agree within a 2 % statistical uncertainty. As indicated 
by the hard band stacked residuals (Fig. |6), EPIC-pn fluxes are 
on average lower than those obtained with XISO, XIS 1 and XIS3 



instruments, with the difference amounting to ~ 10 % (32.5a - ), 
~ 6 % (20.7cr) and ~ 6 % (18.9a") respectively (see Table 0). 
For Suzaku-XIS instruments the object-dependent fraction of the 
PSF scatter of photons from the cluster core to the 3-6 arcmin 
extraction region might affect the observed fluxes. For clusters 
A1795 and A3 112 with the brightest cool cores in our sample, 
we estimated in Sec tion l5"31 that in maximum ~ 15 % of the pho- 
tons in the extraction region may be scattered from the cool core. 
The fraction is expected to be lower for other clusters in our sam- 
ple. Thus, XIS fluxes might be biased high and a correction for 
the PSF effect would thus yield the fluxes more consistent with 
EPIC-pn. 

Since the Chandra-ACIS instrument has been reported to 
yield ~ 10 % higher fl uxes than EPIC-pn in the hard band 
(Nevalain en et alj 12010). we find perfect agreement between 
XISO and ACIS fluxes and XIS1 and XIS3 are in closer agree- 
ment with ACIS than pn. However, correcting XIS fluxes for the 
maximum PSF scatter of ~ 15 % would render all XIS instru- 
ments to closer agreement with pn than ACIS. Thus the com- 
parison if XIS hard band fluxes are consistent with pn or ACIS 
is inconclusive. We also note that hard band XIS / pn cross- 
calibration uncertainties (Fig. |6) yield a maximum systematic 
uncertainty of ~ 10 % to the derived fluxes. 

9. Suzaku OBF contamination 

We examined here whether possible uncertainties in the mod- 
elling of the Suzaku OBF contaminant could explain the prob- 
lems in the cross-calibration of the effective areas of different 
XIS instruments we reported above. 

It has long been observed that the XIS optical blocking fil- 
ter (OBF) suffers from contamination on the spacecraft sidefj 
The contaminant is unknown, but is thought to consist of oxy- 
gen, hydrogen and carbon, and possibly arises from a DEHP- 
like material that has been created by evaporation from the 
spacecraft components. The instrumental and temporal variabil- 
ity of the contaminant in the public Suzaku calibration (which 
used contamination file ae_xi*_contami_20091201.fits) is mod- 
elled using observations of 1E0102-72.3 and RXJ 1856.5-3754. 
The hydrogen-to-carbon number density ratio is assumed to be a 
constant 157.61, whereas the oxygen-to-carbon ratio varies with 
time aJTC| 

^ = 0.27602 x (1 - exp(-(MJD - 53595.34)/341.65)), (2) 
Nc 

where MJD is the Mean Julian Day of the observation. The spa- 
tial variability of the contaminant is measured with monthly in- 
tegrated observations of the sun-lit Earth. The contaminant is 
assumed to have the same composition in all parts of the detec- 
tor, but its column density varies as a function of radius from the 
centre of the detector as 



N c (f,r) 



N c (r,r = 0) 
1 + (r/A(f)) B(,) ' 



(3) 



where r is the radius in arcmin, Nc (f, r — 0) is the column den- 
sity of the contaminant at time t in the centre of the detector and 
Mi) and B(f) are time dependent constants. 

Since the absorption of the contaminant is strongest in the 
soft energy band (see Fig. 1161. possible uncertainties of the col- 
umn densities of the oxygen and carbon in the current public 
calibration of the contaminant may contribute to the soft band 
uncertainties we reported above (Section [72). 



9 See e.g. http://heasarc.nasa.gov/docs/suzaku/analysis/abc/ 

10 See e.g. http://space.mit.edu/XIS/monitor/ 
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Table 4. Fluxes in the hard energy band given by the unabsorbed MEKAL model. 



Hard energy band (2.0 - 7.0 keV) 



Name 


XISO 


XIS1 


XIS3 


pn 




flux" 


flux" 


flux" 


flux" 


A 1060 


1.35 [1.33-1.36] 


1.38 [1.36-1.40] 


1.34 [1.32-1.35] 


1.19 [1.16-1.12] 


A 1795 


0.82 [0.81-0.83] 


0.67 [0.66-0.68] 


0.73 [0.72-0.74] 


1.34 [1.34-1.36] 


A262 


0.63 [0.62-0.64] 


0.61 [0.60-0.62] 


0.60 [0.59-0.61] 


0.55 [0.53-0.56] 


A3 112 


0.46 [0.45-0.47] 


0.45 [0.44-0.46] 


0.46 [0.45-0.47] 


0.43 [0.42-0.43] 


A496 


0.40 [0.40-0.41] 


0.26 [0.25-0.26] 


0.32 [0.31-0.32] 


1.47 [1.45-1.48] 


AWM7 


2.52 [2.49-2.55] 


2.50 [2.47-2.52] 


2.49 [2.46-2.52] 


2.18 [2.17-2.18] 


Centaurus 


2.12 [2.09-2.13] 


1.69 [1.67-1.71] 


1.77 [1.75-1.78] 


2.09 [2.08-2.09] 


Coma 


5.07 [5.06-5.09] 


5.12 [5.10-5.13] 


5.11 [5.10-5.13] 


4.22 [4.20-4.24] 


Ophiuchus 


14.19 [14.16-14.23] 


14.14 [14.11-14.18] 


13.91 [13.88-13.95] 


11.53 [11.50-11.59] 


Triangulum 


4.00 [3.97-4.01] 


3.88 [3.85-3.89] 


3.93 [3.92-3.96] 


3.45 [3.40-3.51] 



Notes. <D) The flux is given in units of 10 13 ergs/cm 2 /s/arcmin 2 and corresponds to the mean surface brightness in the 3-6 arcmin extraction 
region. 



c 
c 




0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 

Energy [keV] Energy [keV] 



Fig. 16. The total (black) transparency of the OBF contaminant 
as a function of energy with contributions of individual elements 
O (blue), C (green) and H (red) for a contaminant with No = 
0.80 x 10 18 crrT 2 and a composition corresponding to an obser- 
vation date of 2006-10-27 (approximately coincident with the 
the mean observation date of our soft band sample), i.e. oxygen- 
to-carbon ratio of 0.20 and hydrogen-to-carbon ratio of 157.61. 



Indeed, if we adopt an estimate for the systematic uncertainty 
from on-axis calibration measurements^ and vary the amount of 
the 2006-10-27 contaminant (see Fig. [TBI) by 1 - 3<x, i.e. vary the 
oxygen column density by 0.5 - 1.5 xlO 17 ctrr 2 while keeping 
the relative abundance of other elements constant, the effective 
area changes increasingly with lower energy (see Fig. [17}. This 
yields the same level of difference (~ 20 % at 0.5 keV) as indi- 
cated by the Suzaku XISO / XIS1 and XIS3 / XIS1 cluster ob- 
servations (see Fig. [TOl . As the change of the oxygen column 
density does not cause any effect in the hard band, agreement 
of the XIS instruments with EPIC-pn in the hard band is main- 



11 from http://space.mit.edu/XIS/monitor/contam/ As we use off-axis 
extraction regions, the uncertainties of the spatial dependence of the 
contaminant will somewhat increase the actual systematic uncertainty 
for our extraction region. The systematic uncertainties of the calibration 
measurements dominate over the statistical uncertainties. 



Fig. 17. The change of the transparency of the 2006-10-27 OBF 
contaminant when the oxygen column density is varied by + 
l.Ocr (dotted lines), 2. Oct (dashed and dotted lines) and 3.0cr 
(dashed lines) while keeping the relative abundance of the other 
elements fixed. 



tained (Section [7. II ). Thus, a possible uncertainty in the oxygen 
and carbon column densities is a feasible cause for the problems 
in the cross-calibration of XIS soft band effective areas. 

9.1. Spectral analysis of OBF column densities 
9.1.1. Method 

We proceeded by testing quantitatively the hypothesis that the 
reported systematic uncertainties of the oxygen column density 
(Section [9]) are responsible for the soft band XIS calibration 
problems (Fig. [TOl . The test includes modifying the XIS effec- 
tive areas with a physically motivated absorption model for the 
contaminant and requiring that all instruments produce an equal 
emission model for a given cluster. 

We adopted the model used by the Suzaku calibration team 
as the absorption model, with absorption cross-sections from 
iBalucinska-Church & McCammonl dl992l) . We determined the 
oxygen-to-carbon ratio for each cluster according to the ob- 
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servation date using the time dependence (Eq. [2]i and used the 
hydrogen-to-carbon ratio of the contaminant in the public cali- 
bration. Q By allowing the oxygen column density to vary, we 
thus varied the total amount of the contaminant, while keeping 
the composition of the model consistent with the contaminant 
implemented in the public calibration. 

We multiplied the MEKAL emission model with the above 
absorption model, while forcing the temperature and metal abun- 
dance equal in all instruments for a given cluster. We allowed the 
normalisation of the emission model for a given cluster to be an 
independent parameter for the different instruments. We created 
auxiliary response files assuming no contaminant, and used these 
to fit the product of the MEKAL model and absorption model. 
While maintaining the emission models equal between all instru- 
ments for a given cluster, the fit to data determines the required 
oxygen column density, i.e. we fit the effective areas simultane- 
ously with the cluster emission. 

9.1 .2. Quality of the modelling 

By regrouping the spectra to uniform binning of 127.75 eV (see 
Section l7.2.1l ). we found that the above modification to XIS ef- 
fective areas does in all cases improve the fits over fits using 
the public calibration (Figs. [81 and IT8l Tables l6l and lA.U . The 
improvement is clearly significant, as is evident from visual in- 
spection of the residuals. However, the standard F-test evaluation 
of the improvement in terms of the decrease of x 1 relative to the 
decrease of the degrees of freedom is not applicable here. While 
the^f 2 in the modified ARF fits decrease, there are more degrees 
of freedom in the modified ARF fits than in the public ARF fits, 
as the number of free parameters in the public ARF fits is 9 (3 
temperatures, 3 abundances and 3 normalisations) while in the 
modified ARF fits the number of free parameters is 8 ( 1 temper- 
ature, 1 abundance, 3 normalisations and 3 values of No).[_J 

9.1 .3. Modification to the contaminant 



XIS joint fit residuals 



When we apply the stacked residuals method (Section [7.2.2t to 
the results with the modified effective areas, we find that the me- 
dian of the ratios R01 and R31 become close to unity (Fig. ITOt. 
Thus these measures of the effective area cross-calibration un- 
certainties improve over those obtained with the public calibra- 
tion. The stacked residuals with modified ARFs are consistent 
with unity in all channels except 0.5 - 0.6 keV for the XIS3 / 
XIS1 ratio. This implies that while the current public H-C-O 
composition with modified column density yields a significant 
improvement over the column density implemented in the cur- 
rent public calibration, there may still be room for adjusting the 
composition of the contaminant and other constituents of the ef- 
fective area. 

In order to evaluate the cluster-based modification to the col- 
umn density of the contaminant implemented in the current cal- 
ibration, we need to work out the values of the oxygen column 
density in our 3-6 arcmin extraction regions in the public cal- 
ibration, since these are not available in the auxiliary response 
files produced by the standard software. 

Currently the radial behaviour of the contaminant (Eq. |3}, 
essential for extended sources like clusters, is only measured 



12 Since the data do not enable the detection of both the C and O edge, 
we are not able to test the calibration of the O-to-C ratio 

13 as spectral fit parameters relating to the sky background (see 
Section l5"Tt are forced equal in all instruments for modified ARF fits, 
the difference in degrees of freedom increase further 




0.5 



1.0 1.5 
Energy [keV] 



2.0 



Fig. 18. The median and the median absolute deviation of the 
residuals (data / best-fit model) of the soft band cluster sample 
for XIS0 (black), XIS1 (red) and XIS3 (blue). The model is an 
absorbed single temperature MEKAL whose temperatures and 
metal abundances are forced equal in all instruments for a given 
cluster, fitted to 0.5 - 2.0 keV band using the effective areas mod- 
ified with the contaminant model. 



for the XIS1 instrument and assumed equal for XIS0 and XIS3. 
This is a possible source of remaining calibration uncertainties 
which we cannot examine with the current data. We divided our 
3-6 arcmin extraction regions into concentric annuli with a 
width of 0.5 arcmin and calculated the No of the public cali- 
bration at the centres of these annuli. We weighted these values 
with th e flux in each annulus u sing the surface brightness pro- 
files oflNevalainenetal.1 d2010h for A 1795, A262 and A3 112; 
iMarkevitch et alj d 19991) for A496 and iTamura et all d2000b for 
A 1060, in order to recover the No values in the current calibra- 
tion. The flux in the annuli remains rather constant with radius 
due to the opposing effects of radially decreasing surface bright- 
ness and radially increasing area of the annuli. Thus, the flux- 
weighted average oxygen column densities are within + 1 % of 
that at a distance of 4.5 arcmin from the centre of the FOV and 
we used this to evaluate the oxygen column densities of the con- 
taminant for our clusters in the public calibration. 

Most of the modified oxygen column densities (see Fig[T9l 
and Table [5]) are within 3cr of the reported uncertainties of the 
values implemented in the public calibration (see Section|9]l. In 
more detail, the cluster analysis indicated that the No values of 
XIS0 are very close to public calibration values up to May 2008 
while XIS1 and XIS 3 require additional oxygen column density 
of the contaminant up to 2 xlO 17 cirT 2 . This implies that the 
soft band effective area shape of XIS0 is currently accurately 
calibrated and that the effective area of XIS1 and XIS3 is over- 
estimated increasingly towards the lowest energies in the public 
calibration for observations up to to May 2008. 

All instruments require a significant increase of the oxygen 
column density for A496 observation on Aug 2008. This effect is 
not evident in the monitoring data of EO1O2-72.0 This may be 
affected by the fact that the Galactic Nh varies by~ 2xl0 20 cm 



within 1° from the cluster centre (Kalberla et al. 



2005b. Also 



allowi ng the column density to be a free parameter, Tanak a et alj 
d2006l) derived a Nh profile for A496 which increases towards 
the cluster centre, reaching a maximum of ~ 6 x 10 20 cm -2 in 



http://www.astro.isas.ac.jP/suzaku/doc/suzaku_td/nodelO.html#contami 
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Table 6. The temperatures, metal abundances and \ 2 Ido f of the 
0.5-2.0 keV band spectral fits when forcing the temperatures 
(kT) and abundances (abund) equal in all instruments for a given 
cluster and allowing the oxygen column density (No) of the con- 
taminant to be a free parameter for all sources and instruments. 
The uncertainties are the statistical ones given at lcr level. 



Name 


kT 


abund 


X z /dof 




[keV] 


Solar 




A 1060 


2.8 [2.7-2.9] 


0.37 [0.35-0.39] 


748.18/613 


A1795 


4.2 [3.0-4.8] 


0.16 [0.12-0.24] 


242.37/246 


A262 


2.1 [2.1-2.2] 


0.55 [0.52-0.59] 


444.91/380 


A3 112 


3.9 [3.7-4.1] 


0.33 [0.29-0.37] 


334.84/326 


A496 


3.4 [3.2-3.6] 


0.29 [0.26-0.36] 


506.83/501 



the 3-6 arcmin annulus, while we used the LAB average of ~ 
4 x 10 20 ctrT 2 in our fits. From consideration of the transmission 
curves (see Fig. [Tol l, we conclude that this discrepancy in the 
Galactic N H can account for the extra contaminant required in 
the fits to A496. 



9.1 .4. New cluster emission models 

Since the required modification of the XIS0 effective area yields 
only minor changes for the oxygen column density of the con- 
taminant (see Fig. [T9T l. the corresponding changes on the XIS0 
temperatures are small (see Fig. [20] and Tables l6l and IA.1) . The 
required higher No for XIS 1 andXIS3 increases the absorption 
and in order to produce the same model prediction, the emission 
models become softer. On average the XIS1 and XIS 3 tempera- 
tures decrease by ~ 27 % (1 1.5cr) and ~ 22 % (7.4<x) respectively 
when compared to the fit with the public calibration. 

Since the soft band temperatures of XIS1 and XIS3 de- 
crease significantly with the modification of the contaminant, 
the Suzaku-XIS v.s. EPIC-pn soft band temperature comparison 
will change from that using the public calibration, presented in 
Section lTXTI The joint XIS0 + XIS1 + XIS3 fit to the soft band 
with the free contaminant yielded on average ~ 12 % (6.2(f) 
lower temperatures compared to EPIC-pn soft bands tempera- 
tures (see Fig. |2TT>. As noted in Section 1331 the PSF scattering 
from cool cores may bias XIS best-fit temperatures low by 7 
% in the worst case; however, we expect this bias to be much 
smaller in most cases. If corrected for the PSF scatter, the XIS 
soft band temperatures might increase by a few %, i.e., become 
closer to the EPIC-pn values. Since the Chandra-ACIS instru- 
ment has been found to yi eld ~18% higher tempe ratures than 
EPIC-pn in the soft band dNevarainen et alJl2010l) . our results 
indicate that the XIS instruments yield a better soft band tem- 
perature agreement with the EPIC-pn instrument than with the 
Chandra-ACIS instrument. The data exhibit quite a lot of scatter: 
three of the clusters yield consistent values between Suzaku-XIS 
and EPIC-pn, while the temperatures of two of the clusters differ 
by 4.9cr and 7. lcr respectively. 

With the introduction of the modified contaminate the 
Suzaku-XIS 1,3 / EPIC-pn stacked soft band residuals at the 
lowest energies become consistent with unity, while XIS0 and 
pn remain consistent in both cases (compare Figs. ITTI and |22Ti. 
However, the scatter is larger (10-20 %) when using the mod- 
ified contaminate, and thus the remaining energy-dependent un- 
certainties, which amount to 12 % temperature differences (see 
above), are not significantly detected with the current data. 



A 
I — 

V 




2 4 6 

XIS joint fit soft bond T [keV] 

Fig. 21. Upper panel: The best-fit soft band temperatures (as- 
terisks) and lcr uncertainties (boxes) for EPIC-pn and Suzaku 
XIS instruments using the modified XIS effective area. The 
solid black line is an identity line, drawn as a reference. Lower 
panel: The relative temperature differences fj (diamonds) be- 
tween EPIC-pn and Suzaku, and their lcr uncertainties. The dot- 
ted and dashed lines show the weighted means of the relative 
temperature difference and corresponding uncertainties. 
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Fig. 22. The median + median absolute deviation of the soft en- 
ergy band stacked residuals of XIS0 / pn R m o pn (top panel), XIS 1 
/ pn R m i pn (middle panel) and XIS 3 / pn R m 3 pn (bottom panel) 
using the contamination-modified effectve areas for XIS instru- 
ments. The data are regrouped to uniform binning of 127.75 keV. 
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2006 



2007 



2008 



Fig. 19. The upper panel shows the time de- 
pendence of the oxygen column densities and 
1.0 cr statistical uncertainties of the contami- 
nant fitted to the cluster data (diamonds) and 
the ones implemented in the public calibration 
at 4.5 arcmin distance from the centre of the 
FOV (lines) for XIS0 (black), XIS1 (red) and 
XIS3 (blue). The lower panel shows the time 
dependence of the difference of these. 



Table 5. The oxygen column density of the 0.5-2.0 keV band spectral fits when forcing the temperatures (kT) and metal abundances 
(abund) equal in all instruments for a given cluster and allowing the oxygen column density (No) of the contaminant to be a free 
parameter for all sources and instruments. A No describes the modification to the oxygen column density of the contaminant. The 
uncertainties are the statistical ones given at lcr level. 







XIS0 




XIS1 






XIS 3 




Name 


N 


A N 




N 


A N 




N 


A N 






[10 18 cm" 2 ] 


[10 18 cm 


- 2 ] 


[10 18 cm" 2 ] 


[10 18 cm 


- 2 ] 


[10 18 cm" 2 ] 


[10 18 cm 


- 2 ] 


A 1060 


0.14 [0.12-0.15] 


0.04 [0.02 - 


0.04] 


0.20 [0.19-0.21] 


0.10 [0.09- 


0.10] 


0.37 [0.36-0.38] 


0.11 [0.10- 


0.12] 


A 1795 


0.15 [0.13-0.17] 


0.02 [0.00 - 


0.04] 


0.24 [0.23-0.26] 


0.10 [0.08- 


0.11] 


0.44 [0.41-0.47] 


0.11 [0.08- 


0.14] 


A262 


0.85 [0.81-0.89] 


0.00 [-0.05 - 


0.03] 


1.09 [1.07-1.12] 


0.20 [0.18- 


0.23] 


1.31 [1.26-1.36] 


0.16 [0.11- 


0.20] 


A3112 


1.06 [0.97-1.10] 


-0.06 [-0.15- 


-0.02] 


1.08 [1.05-1.13] 


0.10 [0.07- 


0.15] 


1.20 [1.15-1.29] 


-0.02 [-0.07 - 


- 0.07] 


A496 


1.33 [1.28-1.38] 


0.16 [0.10- 


0.20] 


1.30 [1.26-1.33] 


0.30 [0.26 - 


0.33] 


1.50 [1.44-1.55] 


0.28 [0.22 - 


0.33] 



Table 7. The relative difference (//) and significance (sig.) of the 
temperatures and fluxes obtained with different instrument com- 
binations. 



Instruments 


1 hard 


1 soft 


flux hard " 1 






sig. d 


» c 


sig. * 


ft c 


sig. 


ACIS/pn c 


-l 


0.6 


18 


8.6 


11 


24.7 


pn/XISO 


5 


4.2 


14 


6.2 


-5 


1.1 


pn/XISl 


2 


1.7 


-18 


7.4 


-1 


0.2 


pn/XIS3 


6 


4.9 


-9 


3.0 


-1 


0.2 


pn/XISO 13 f 






12 


6.2 






XIS0/XIS1 


-5 


5.6 


-29 


11.2 


4 


1.0 


XIS0/XIS3 


-1 


0.7 


-23 


7.4 


4 


1.0 


XIS 1 / XIS 3 


5 


4.9 


9 


2.9 





0.0 



Notes. The values correspond to the best-fit single temperature models 
in (o) 2.0 - 7.0 keV and {b) 0.5 - 2.0 keV energy bands using the public 
calibration, i.e. for Suzaku this refers XIS CALDB version 20110608 
and XRT CALDB version 20080709 or using the modification to the 
XIS contaminant <0 as explained in the text. (c) fi gives the weighted 
mean of the relative difference between the measured temperatures in 
percentages; (rf) sig. gives the stat istical significance of th e temperature 
difference in terms of cr ; (e) from Nevalainen et al. (20lo|). 



10. Conclusions 

We used a sample of galaxy clusters to perform a study on the 
accuracy of the cross-calibration of the energy dependence and 
the normalisation of the effective area of XIS0, XIS1 and XIS3 
instruments onboard Suzaku and EPIC-pn instrument onboard 
XMM-Newton in 3 - 6 arcmin annular extraction regions. The 
Suzaku results apply to the XIS CALDB versions 201 10608 (i.e. 
public calibration at the writing of this paper) and 20110201, 
and to XRT CALDB version 20080709 for observations carried 
out in the time period of years 2005 - 2008. The XMM-Newton 
results refer to SAS version xmmsas_201 10223_1801-1 1.0.0 and 
the latest calibration information in June 2011. The results are 
summarised in Fig. [23] and Table [7] 

As the shape of the spectrum determining the temperature of 
the cluster depends on the energy dependence effective area, and 
the flux depends on the normalisation of the effective area we 
performed the study by comparing the temperatures, fluxes and 
best-fit emission models of the clusters measured with different 
instruments in the 0.5 - 2.0 and 2.0 - 7.0 keV bands. 

Comparison of the Suzaku-XIS instruments with each other 
yielded that XIS0 and XIS3 produced rather consistent hard 
band temperatures which are systematically lower than those ob- 
tained with XIS 1 by ~ 5 %. The Suzaku-XIS hard band tempera- 
tures are systematically lower than those obtained with EPIC-pn 
instrument of XMM-Newton, but only by 2 - 6 % in average. 
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Fig. 20. Upper panel: The best-fit soft band temperatures (asterisks) and 1 cr uncertainties (boxes) for Suzaku XIS instruments using 
the modified XIS effective area (joint fit) and public calibration. The solid black line is an identity line, drawn as a reference. Lower 
panel: The relative temperature differences fj (diamonds) between the Suzaku XIS instruments, and their lcr uncertainties. The 
dotted and dashed lines show the weighted means of the relative temperature difference and corresponding uncertainties. 
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Fig. 23. The average relative difference (dia- 
monds) ± the error of the mean of the temper- 
atures (solid line) and fluxes (dotted line) us- 
ing the public calibration in the soft band (left 
side of the plot) and in the hard band (right 
side of the plot). Comparison of the pn and XIS 
soft band temperature using the modification to 
the contaminant is shown with diamonds and 
dashed line. 



Some fraction of this difference could be caused by Suzaku's 
PSF scattering from the cool core, which should be less than ~ 
6 % in the worst case. Since EPIC-pn hard band temperatures 
are consistent with Chandra-ACIS and BeppoSAX-MECS tem- 
peratures, and EPIC-pn bremsstrahlung temperatures agree wit h 
Fe XXV/XXVI line ratio temperatures dNevalainen et al.ll2010l) . 
we suggest that the energy dependence of XIS instruments is 
quite accurately calibrated in the hard band. The remaining cal- 
ibration uncertainties may yield effects on the XIS hard band 
temperatures at ~ 6 % level. 

In the soft 0.5 - 2.0 keV energy band the XIS1 and XIS3 in- 
struments yielded consistent temperatures when using the public 



calibration, while XISO yielded temperatures which are lower by 
an average of 29 % ( 1 1 .2<x) and 23 % (7.4<x) from those obtained 
with XIS1 and XIS 3, respectively. Comparison of the residu- 
als showed that XISO effective area is underestimated or XIS 1 
and XIS3 effective area is overestimated towards lower energies, 
with the difference increasing to ~ 20 % at 0.5 keV. 

We tested the assumption that the above discrepancies in the 
effective area in the soft band are due to uncertainties of the im- 
plemented column density of the optical blocking filter (OBF) 
contaminant. We achieved this by allowing the OBF contami- 
nant to have a column density different from that in the public 
calibration and fitting jointly the XISO + XIS1 + XIS3 soft band 
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spectral data, while forcing the temperatures and metal abun- 
dances equal in all instruments for a given cluster. The resulting 
fits to the data were significantly better than the ones using the 
public calibration. Assuming the composition of the contami- 
nant as in the public calibration, we found that the cluster data 
required an addition to the column density of the contaminant of 
XIS1 and XIS3 amounting to amaximumANo of 2 xlO 17 cm" 2 , 
i.e. at a 4<x level of the reported uncertainties of the No measure- 
ments in the public calibration. The column density of the XISO 
contaminant implemented in the public calibration is consistent 
with that obtained in our spectral analysis of clusters. Thus, the 
analysis implies that the effective area of XISO in the soft band 
is more accurately modelled in the public calibration than that of 
XIS1 and XIS3, for observations during 2005 - 2008. 

The XIS soft band temperatures obtained with the modified 
effective area are lower by ~ 12 % from those obtained with 
EPIC-pn instrument onboard XMM-Newton. Considering that 
the PSF scatter from the cool cores may bias XIS soft band 
temperatures low, the PSF-corrected XIS soft band tempera- 
tures would become in better agreement with those obtained 
using EPIC-pn. However, Chandra-ACIS instrument yields ~ 
1 8 % higher temperat ures than EPIC-pn in the soft band 
( dNevalainen et alJ20 10)). Thus, our results indicate that the XIS 
instruments yield a better soft band temperature agreement with 
the EPIC-pn instrument than with the Chandra-ACIS instrument. 

The hard band fluxes derived with Suzaku-XIS instruments 
are higher than those derived with EPIC-pn by ~ 6 - 10 %, 
whereas the Chandra-ACIS instrumen t yields ~ 10 % highe r 
fluxes than EPIC-pn in the hard band (Nev alainen et al.l l2010). 
However, as Suzaku-XIS fluxes might be biased high by a maxi- 
mum of ~ 15 % due to PSF scatter and Suzaku-XIS and EPIC-pn 
hard band data shows that the derived fluxes have a maximum 
systematic uncertainty of ~ 10 %, comparison of Suzaku-XIS to 
EPIC-pn and Chandra-ACIS fluxes proved inconclusive. 
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Appendix A: Spectral fits 

Here we show the data and the best-fit single temperature MEKAL models for 
EPIC-pn (Fig.|AT]l, XISO (Fig.|A2]l, XIS1 (Fig.[0 and XIS3 (Fig.fA4) using 
the public calibration. The spectral parameters are listed in Table lA. il 
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Table A.l. Best-fit parameters of the Suzaku XIS and XMM-Newton EPIC-pn spectral fits using the public calibration 







XIS0 






XIS1 




Name 


kT 


Abund 


X 2 ldof 


kT 


Abund 


X 2 /dof 




[keV] 


[Solar] 




[keV] 


[Solar] 




Soft energy band (0.5 - 2.0 keV) 


A 1060 


2.8 [2.7-2.9] 


0.43 [0.39-0.48] 


229.51/170 


4.7 [4.6-4.9] 


0.98 [0.91-1.05] 


413.27/275 


A 1795 


4.2 [3.8-4.6] 


0.26 [0.15-0.35] 


73.02/66 


11.0 [9.7-13.5] 


0.03 [0.00-0.32] 


125.64/108 


A262 


2.0 [2.0-2.1] 


0.51 [0.47-0.62] 


166.66/108 


2.4 [2.3-2.5] 


0.85 [0.80-0.91] 


215.40/159 


A3 112 


3.2 [3.0-3.6] 


0.21 [0.14-0.29] 


84.00/86 


4.7 [4.3-5.0] 


0.61 [0.49-0.69] 


165.84/140 


A496 


4.0 [3.8-4.2] 


0.45 [0.37-0.55] 


170.79/146 


5.5 [5.3-5.8] 


0.85 [0.75-0.95] 


289.28/205 


Hard energy band (2.0 - 7.0 keV) 


A 1060 


3.5 [3.4-3.6] 


0.44 [0.40-0.48] 


109.55/106 


3.4 [3.3-3.5] 


0.53 [0.48-0.58] 


135.32/121 


A 1795 


6.6 [6.2-7.0] 


0.40 [0.34-0.47] 


43.04/41 


5.9 [5.5-6.2] 


0.39 [0.33-0.46] 


56.73/51 


A262 


2.5 [2.5-2.6] 


0.45 [0.40-0.52] 


43.91/50 


2.5 [2.4-2.5] 


0.52 [0.45-0.60] 


76.81/57 


A3 112 


5.4 [4.7-5.7] 


0.45 [0.40-0.51] 


78.00/60 


5.5 [5.3-5.9] 


0.40 [0.35-0.46] 


69.98/74 


A496 


4.2 [4.1^4.3] 


0.49 [0.45-0.53] 


149.27/119 


4.4 [4.3-4.6] 


0.51 [0.47-0.56] 


130.02/129 


AWM7 


3.9 [3.8-1.0] 


0.58 [0.53-0.63] 


118.18/97 


4.0 [3.9-4.1] 


0.71 [0.66-0.77] 


112.21/105 


Centaurus 


3.6 [3.6-3.7] 


0.81 [0.77-0.86] 


171.03/159 


3.8 [3.7-3.9] 


0.82 [0.77-0.87] 


223.60/177 


Coma 


8.3 [8.2-8.4] 


0.36 [0.35-0.38] 


1152.77/1124 


9.0 [8.8-9.1] 


0.37 [0.36-0.38] 


1155.40/1104 


Ophiuchus 


9.9 [9.8-10.0] 


0.47 [0.46-0.49] 


1337.33/1242 


10.7 [10.5-10.8] 


0.46 [0.44-0.47] 


1282.61/1187 


Triangulum 


9.5 [9.3-9.8] 


0.31 [0.29-0.33] 


512.69/537 


10.3 [10.0-10.6] 


0.35 [0.32-0.38] 


547.64/536 






XIS3 






pn 




Name 


kT 


Abund 


X 2 ldof 


kT 


Abund 


X 2 /dof 




[keV] 


[Solar] 




[keV] 


[Solar] 




Soft energy band (0.5 - 2.0 keV) 


A 1060 


3.8 [3.5-3.9] 


0.76 [0.68-0.83] 


239.18/159 


2.9 [2.8-2.9] 


0.48 [0.45-0.50] 


327.00/290 


A 1795 


6.5 [5.6-7.5] 


0.51 [0.28-0.82] 


61.45/61 


4.2 [4.1-4.3] 


0.20 [0.18-0.23] 


331.04/296 


A262 


2.5 [2.4-2.6] 


0.75 [0.67-0.83] 


102.89/102 


2.5 [2.4-2.5] 


0.71 [0.67-0.78] 


313.94/240 


A3 112 


3.9 [3.7-4.5] 


0.32 [0.24-0.45] 


90.14/89 


3.7 [3.6-3.8] 


0.24 [0.22-0.28] 


353.08/296 


A496 


4.9 [4.6-5.1] 


0.65 [0.53-0.76] 


131.56/141 


4.9 [4.8-5.0] 


0.72 [0.69-0.78] 


472.49/297 


Hard energy band (2.0 - 7.0 keV) 


A1060 


3.4 [3.3-3.4] 


049 [0.44-0.53] 


108.91/104 


3.5 [3.3-3.6] 


0.42 [0.37-0.47] 


124.03/167 


A1795 


6.6 [6.2-7.1] 


0.34 [0.28-0.40] 


50.94/42 


7.3 [7.1-7.5] 


0.35 [0.32-0.38] 


324.32/400 


A262 


2.4 [2.3-2.4] 


0.41 [0.36-0.48] 


46.48/48 


2.5 [2.4-2.6] 


0.50 [0.41-0.60] 


38.66/78 


A3 112 


5.3 [5.0-5.5] 


0.45 [0.39-0.50] 


77.74/63 


5.9 [5.6-6.2] 


0.42 [0.38-0.46] 


311.65/334 


A496 


4.6 [4.2-4.8] 


0.43 [0.39-0.47] 


122.60/114 


5.0 [4.9-5.1] 


0.46 [0.43-0.49] 


277.85/353 


AWM7 


4.1 [3.9-4.2] 


0.59 [0.55-0.64] 


109.10/96 


4.0 [4.0-4.0] 


0.57 [0.56-0.59] 


883.38/806 


Centaurus 


3.6 [3.5-3.7] 


0.81 [0.77-0.86] 


176.12/157 


3.8 [3.8-3.9] 


0.71 [0.69-0.73] 


814.58/655 


Coma 


8.6 [8.5-8.7] 


0.35 [0.34-0.36] 


1356.87/1151 


9.0 [8.8-9.3] 


0.35 [0.33-0.38] 


505.40/531 


Ophiuchus 


9.9 [9.7-10.0] 


0.47 [0.46-0.48] 


1392.54/1257 


10.4 [10.2-10.7] 


0.41 [0.38-0.43] 


677.17/703 


Triangulum 


9.8 [9.6-10.1] 


0.32 [0.30-0.35] 


573.44/524 


11.3 [10.3-12.7] 


0.40 [0.29-0.50] 


62.88/106 
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Fig. A.l. The soft band (upper panel) and hard band (lower panel) XMM-Newton EPIC / pn spectra (crosses) and best-fit single 
temperature independent fits (solid lines) for the cluster sample. The normalisations of the spectra are adjusted for plot clarity and 
do not reflect the relative brightness of the clusters. 
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Fig. A.2. The soft band (upper panel) and hard band (lower panel) Suzaku XISO spectra (crosses) and best-fit single temperature 
independent fits (solid lines) for the cluster sample. The normalisations of the spectra are adjusted for plot clarity and do not reflect 
the relative brightness of the clusters. 
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Fig. A.3. The soft band (upper panel) and hard band (lower panel) Suzaku XIS 1 spectra (crosses) and best-fit single temperature 
independent fits (solid lines) for the cluster sample. The normalisations of the spectra are adjusted for plot clarity and do not reflect 
the relative brightness of the clusters. 
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Fig. A.4. The soft band (upper panel) and hard band (lower panel) Suzaku XIS3 spectra (crosses) and best-fit single temperature 
independent fits (solid lines) for the cluster sample. The normalisations of the spectra are adjusted for plot clarity and do not reflect 
the relative brightness of the clusters. 
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